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Abstract 



We discuss implications of improved upper bounds on the |AL| = 2 pro- 
cesses (i) K + — > n~fj, + fj, + , from an experiment at BNL, and (ii) pT — > e + 
conversion, from an experiment at PSI. In particular, we address the issue of 
constraints on neutrino masses and mixing, and on supersymmetric models 
with i?-parity violation. 
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At present there are increasingly strong indications for neutrino oscillations and hence 
neutrino masses and lepton mixing from the solar neutrino deficiency and atmospheric neu- 
trino anomaly [l] . The existence of lepton mixing means that lepton family number is not a 
good symmetry. Majorana neutrino masses occur generically, and violate total lepton num- 
ber L by |AL| = 2 units. However, so far, in contrast to the data suggesting lepton mixing, 
experimental searches for the violation of total lepton number have only set limits. Among 
these are searches for the |AL| = 2 processes (i) neutrinoless double beta (Oz/2/?) decay of 
nuclei and (ii) pT — ► e + conversion in the field of a nucleus. A third class of |AL| = 2 
processes includes the decays K + — > where = e + e + , fi + e + , or In 

a previous work we considered these decays and, from a retroactive data analysis, set the 
first upper limit on one of them, namely [|]|| 

BR(K+ -> 7v-/2 + i2 + ) < 1.5 x 1(T 4 (90% CL). (1) 

In [|J we also noted that rare K decay experiments at BNL could greatly improve this limit 
and proposed a search for K + — > 7i~fi + fi + 0. Among these experiments was BNL E865, 
which was searching for, and has now set a stringent upper limit on, the decay K + — > 7i + fi + e~ 
0. This experiment has also recently obtained the 90 % CL upper limit 

BR(K + -> 7T-fi + fi + ) < 3.0 x 10~ 9 . (2) 

In the present paper we discuss the implications of this limit. In the context of current 
bounds on neutrinoless double beta decay and fi~ — > e + conversion, we shall also consider 
the implications of two other 90 % CL limits on |AL| = 2 decays from E865 [Q: 

BR(K + -> 7r"e + e + ) < 6.4 x 10" 10 (3) 

and 

BR(K + -> 7r"/i + e + ) < 5.0 x 10~ 10 . (4) 

We first discuss physics sources for these decays, concentrating on massive Majorana 
neutrinos and i?-parity-violating supersymmetric (SUSY) theories. In a modern theoretical 
context, one generally expects nonzero neutrino mass terms, of both Dirac and Majorana 
type. Let us denote the left-handed flavor vector of SU(2) x U(l) doublet neutrinos as 
Vl = (uejU^u^L and the right-handed vector of electroweak-singlet neutrinos as iV^ = 
(Ni, .., N ns )ji. The Dirac and Majorana neutrino mass terms can then be written compactly 
as 

1 f M L M D \ ( v% \ , 

- C m = -{ VL JVJ ( (M ; )T m r ){ n R r )+ h.c. (5) 

where Ml is the 3x3 left-handed Majorana mass matrix, Mr is a n s x n s right-handed 
Majorana mass matrix, and Mp is the 3-row by n s -column Dirac mass matrix. In general, 
all of these are complex, and (M^) T = Mr, , (M^) T = Mr. The diagonalization of the 
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matrix in eq. (|5]) then yields 3 + n s mass eigenstates, which are generically nondegenerate 
Majorana neutrinos (degeneracies in magnitudes of eigenvalues can yield Dirac neutrinos). 
Writing the charged current in terms of mass eigenstates as J\ = £lJ\Vl, one has, in par- 
ticular, Up = X]j=™ s UfijVj. The seesaw mechanism naturally yields a set of three light 
masses for the three known neutrinos, generically of order m v ~ mj^/MR, and n s very large 
masses generically of order tur, for the electroweak-singlet neutrinos, where mn ~ Mew and 
itir » Mew denote typical elements of the matrices Me and Mr, and Mew — 250 GeV 
is the electroweak symmetry breaking scale. However, if one tries to fit all current neutrino 
experiments, including not only the solar and atmospheric, but also the LSND data, then it 
is necessary to include electroweak-singlet ("sterile") neutrinos with masses << Mr in order 
to achieve an acceptable fit. In this case, the weak eigenstate may contain significant 
components of mass eigenstates beyond the usual three light ones Uj, j = 1,2,3, and, it is 
a priori possible that some of these mass eigenstates might have masses lying in the theo- 
retically "disfavored" intermediate range m,D/m R « m v . « m R (subject to both particle 
physics and astrophysical/cosmological constraints). Whether or not such intermediate-mass 
neutrinos exist is ultimately an empirical question that must be settled by experiment. It 
is therefore of continuing interest to address constraints from data on neutrino masses in 
this intermediate mass region. Since current indications from solar and atmospheric data 
suggest neutrino masses in the seesaw-favored region and since, independent of this, neutrino 
masses in the range of a few to several hundred MeV can suppress large-scale structure for- 
mation and hence may be disfavored by cosmological constraints ||, we concentrate mainly 
on neutrino masses smaller or larger than this range here. 

There are two types of lowest-order graphs involving massive neutrinos that contribute 
to the decay K + — > n~£ + £' + , where £ + £' + = /i + e + , or e + e + , as shown in Fig. 1. 
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FIGURES 



FIG. 1. Graphs involving massive Majorana neutrinos that contribute to K + ir~£ + £' + , where 
£ + ,£' + = H + e + , or e + e + . In the case of identical £ + and £' + , it is understood that the contributions 

are from the diagrams minus the same diagrams with the outgoing antilepton lines crossed. 



For the K + — > 7r decay, the s-channel diagrams yield, in standard notation, 

Amp{2{a)) = 2G 2 F f K f 7r (V ud V us )* ^{U^Yp^jp^ [if (p^jy) - Lfip^pj] (6) 

j 

where 

Lf p,) = m Uj [q 2 - ml^viP.h^PRv'M (7) 

where q denotes the momentum carried by the virtual Uj, and Pr = (1 + 7s)/2 is a right- 
handed chiral projection operator. The i-channel graphs cannot be evaluated so easily, 
because the hadronic matrix element that occurs, 

M^'V'^' 1 (tT I [d L {y) lp u L {y)\ [s L (x) la u L {x)\ \ K + ) (8) 

cannot be directly expressed in terms of measured quantities, unlike the matrix elements 
(0\s L , ~f a u L \K Jr ) and (ir~ \d,Lj/3U L \0) from the first graph. In the limits where m 2 . is much 
smaller or larger than the mag nitude of the typical g 2 , |(g 2 ) a „ e | ~ O((10 2 MeV) 2 ), the 
propagator factor simplifies: 



va v . 



"3 



m v ./{q 2 ) ave iim V] «\(q 2 ) c 
q2 ~ m ^ ) -l/m„. ifm„. » \{q 2 ) c 



(9) 



Since the contribution of each virtual Vj is accompanied by a complex factor (U*j) , it 
is possible for these contributions to add constructively or destructively. Because of the 
possibility of such cancellations, one cannot put an upper limit on neutrino masses or lepton 
mixing matrix coefficients from an upper bound on the decay K + — ► Tr~p + p, + . (A similar 
remark applies to K + — > 7i~fi + e + and K + — > 7r~e + e + since cancellations can also occur 
among contributions of various z/j's to the respective amplitudes for those decays.) The 
same comment is well-known in the case of neutrinoless double beta decay; for example, in 
the light-mass region, the lower limit on half-lives for Qv2f3 transitions places an upper limit 
on ^2jU 2 jm Uj , not on \U e j\ or m Uj themselves. 
With these inputs, we estimated [|J 

BR(K + - 7r> + /i + ) ~ 10-( 13±2 Vj^^(m„./(100 MeV) * (10) 
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where 



z if z « 1 
1/z if z » 1 



(11) 



and the factor r 
K + — > 7r _ e + e" 



^ 0.2 is a relative phase space factor (normalized relative to the decay 
. If z ~ 0(1), one must, of course, retain the exact propagator. In Fig. 



? < q 2 < 



m K + - m. 



i.e., 245 < 



< 388 



1(a), the range of (timelike) q is (m w + + m A 
MeV. Hence, if there exists a neutrino with m u . in this range, 245 < m v . < 388 MeV, then 
q 2 — m 2 ,. can vanish, leading to a resonant enhancement of the amplitude for this graph [ |I0|| . 

In the following, we assume for simplicity that a single mass eigenstate Uj dominates 
the sum in (|1(]) but recall our remark above concerning the possibility of cancellations and 
note that it is straightforward to generalize our discussion to the case where there are several 
comparable contributions. If m v . « rriK, then, using the new upper bound (|2[) and the most 
conservative choice for the numerical prefactor in eq. dig), (i.e., taking 1(H 13±2 ) -> 1(T 15 ), 
we obtain 



\U, 



i-ij i 



m. 



100 MeV 



< 4 x 10' 



BR(K^ 



71 fl 



3 x 10- 9 



1/2 



On the other hand, if m u . >> m^, we obtain 



\U, 



i-ij i 



100 MeV 



< 4 x 10' 



BR{K + -> 7r-/i+/x H 
3 x 10- 9 



1/2 



(12) 



(13) 



Since m v . « m K in order for ( |T2"D to hold, and since |C/ W -| < 1 by unitarity, the bound ( |I~2"D 
does not place a significant restriction on either of these quantities. Similarly, in the heavy 
neutrino mass region, since m u . » tjik in order for the bound (|T3| ) to hold, it does not 
place any restriction on m v . or \U^j\ in this mass region |14]]. This is not to say, however, 



that it is not worthwhile to search further for the decay K + — > ti~ since it constitutes 
a testing ground for violation of total lepton number that is quite different from the usual 
searches for neutrinoless double beta decay of nuclei. 

We also give an update of the limit on the decay K + — > 7r~fi + e + . In [|J we obtained an 
indirect upper limit on this decay by observing that the leptonic part of the amplitude for this 
decay is related by crossing to the leptonic part of the amplitude for the conversion process in 
the field of a nucleus (Z, A): fx~ + (Z, A) — ► e + + {Z — 2, A). At the time of 0], the best upper 
limit on this conversion process was a(fi~+Ti — > e + +Ca)/<j(n~+Ti — > u^+Sc) < 1.7xl0~ 10 

x 10~ 9 . 



from a TRIUMF experiment [15|, and we inferred that BR(K^ 
The current best bound, from a PSI experiment, is 



7T 



< 



few 



cr(/i + Ti — > e H 



° a) < 1.7 x 10- 12 



cr(fi + Ti — > + Sc) 
Using this new bound, we conservatively infer that 

BR(K + -> 7T-fi + e + ) < few x 10" 11 



(90 % CL) . 



(14) 



(15) 



5 



As was noted in this is an indirect limit since it requires a theoretical estimate of the 
hadronic matrix element as input. Our indirect bound flT5| ) is more stringent than the E865 
limit (H), but the latter limit is still valuable since it is direct. 

One can obtain an indirect upper limit on the decay K + — * 7r - e + e + from the existing 
upper limit on neutrinoless double beta decay, which, for light neutrinos, gives | J2j U^m^ | < 
0.4 eV (depending on the input used for the nuclear matrix elements) [[1711 . Using the same 
method as above, we obtain an upper limit many orders of magnitude less than the direct 
limit (|). 

One can also consider AL = 2 decays of D and B mesons. Here only rather modest 
direct upper limits of order 10~ 3 to 10~ 4 have been set on the branching ratios ]l!|. A 



similar comment applies to |AL| = 2 hyperon decays, on which we previously set upper 
limits |18 . 

We next proceed to discuss constraints on i?-parity violating (RPV) SUSY models, im- 
parity may be defined as R = (— i) 3B + L + 2S ' ) where B, L, and S refer to the baryon and lepton 



numbers and to the spin of the particle [Bl . Although i?-parity was originally hypothesized 



in order to prevent intolerably rapid proton decay in SUSY models, one can achieve the 
same end by imposing weaker global symmetries that forbid terms of the form U^D C -D k in 
the superpotential (where the subscripts i, j, k here are generation indices and we follow the 
usual convention of writing the holomorphic operator products in terms of left-handed chiral 
superfields) while still allowing the i?-parity violating terms 

W RPV = X ijk LiLjE c k + X' ijk LiQjDl + . (16) 

These terms violate total lepton number and, in general, also lepton family number. A recent 
review of i?-parity violating SUSY models is ||20|| . The second term in (|1^) yields several 
contributions to K + — > 7r~ / u + / u + , shown in Figs. 2(a,b), where jl and x° denote the scalar 
muon and neutralino. 




6 



8 * ' ' £' + 
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FIG. 2. Graphs that contribute to K + — > 7r~£ + i" + , where = /i + /J + , /i + e + , or e + e + , in supersym- 

metric theories with R p violation. In the case of identical £ + and £ /+ , it is understood that the contributions 
are from the diagrams minus the same diagrams with the outgoing antilepton lines crossed. 

Note that there may be neutrino-neutralino mixing in RPV theories, even if at some 
mass scale one rotates the terms KiLiH u to zero. A third type of diagram is shown in 
Fig. 2(c). As noted, for K + — > tt~u + u + , each diagram is accompanied by minus the same 
diagram with the outgoing u + lines crossed. The ix ^ vertices are cx \J g 1 + g' 2 , where g 
and g' are the SU(2) and U(l)y gauge couplings, while the udp, and sup, vertices are cx A' 211 
and A' 212 , respectively and the ude and sue vertices are cx A' m and A' 211 , respectively. In the 
third type of diagram, there can be a gluino on the internal line, with dgd and u c gu vertices 
proportional to the strong coupling g s ; or there can be a neutralino on the internal line, 
with vertices as given above. Bounds on RPV couplings are model-dependent, but typical 
current upper bounds on A' 211 , A 212 , A' m , and A' 211 are < O(0. 1) [pOfl . Using these inputs, 
we find that these i?-parity violating contributions could be much larger than those from 
massive neutrinos and lepton mixing, but are still expected to be small compared with the 
upper limit @: 



BR(K + -> vr>VW S 10- 16 (A 211 A 212 ) 2 (^^) 10 . (17) 

v m S usY ' 



For the purpose of this rough estimate, we have taken the masses of the various superpartners 
u, d, p, x°, and g to be comparable and denoted this mass scale as msusY ~ Mew- Note that, 
given the lower bounds on the masses of e and jl or order 100 GeV, no resonance is possible in 
the amplitude of Fig. 2(a). One could, of course, take a particular SUSY parameter set and 
perform the calculation for this set; however, there is a large range of variation in possible 
superpartner masses as well as allowed values of other relevant parameters such as tan j3 (as 



illustrated, e.g., by the parameter sets used in plf), so we deliberately keep our estimate 
general. Hence, it appears that the limit @ does not strongly constrain possible RPV 
SUSY theories. The constraints on i?-parity violating models from yT — > e + conversion and 
neutrinoless double beta decay are also more stringent than those from the limits (Q) and 

(0). 

Thus, while neutrinoless double beta decay and \x — > e + conversion are the most sensitive 
ways to search for |AL| = 2 transitions with |AL e | = 2 and AL e = AL M = ±1, respectively, 
the decay K + — > tt~u + u + is, at present, the best way to search for |AL| = 2 transitions 
with |AL M | = 2. It is therefore worthwhile to estimate the potential of future K + decay 
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experiments to probe for this decay to lower values of branching ratio. In particular this 



might be undertaken by the CKM experiment planned at Fermilab ||22|| . The proposal for this 



experiment anticipates a statistical sensitivity of ~ 10~ 12 /event for the decay K + — > 7i + uu. 
While it would require a substantial, and as yet unplanned, effort to design and build a 
trigger to search for K + — > 7i~fi + fi + , it appears |22| that this experiment might be able to 



reach a level near to 10~ 12 in branching ratio and thus improve substantially on the already 
impressive upper limit on BR(K + — > 7i~fi + fi + ) from E865 at Brookhaven. We also suggest 
that the planned experiment MECO at BNL [^3fl , which anticipates searching for — > e~ 
conversion below the level ~ 10~ 16 relative to fi~ capture, should also undertake a search for 
fi~ — > e + conversion. 
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government retains a non-exclusive royalty-free license to publish or reproduce the published 
form of this contribution or to allow others to do so for U.S. government purposes. 
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